Phase Diagrams and Crossover in 
Spatially Anisotropic d=3 Ising, XY Magnetic and Percolation Systems: 
Exact Renormalization-Group Solutions of Hierarchical Models 



o 
o 

(N 
5h . 

< 
in 

(N 



CO 

c3 



I 

C 

O 

o 



> 

t> 

O 

IT) 

o 



c 

o 
o 



X 
S3 



Aykut Erba§, 1 Ash Tuncer, 1 Burcu Yiicesoy, 1 and A. Nihat Berker 1-3 
Department of Physics, Istanbul Technical University, Maslak 34469, Istanbul, Turkey, 
2 Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, U.S.A., and 
3 ' Feza Giirsey Research Institute, TUBITAK - Bosphorus University, Cengelkdy 81220, Istanbul, Turkey 

Hierarchical lattices that constitute spatially anisotropic systems are introduced. These lattices 
provide exact solutions for hierarchical models and, simultaneously, approximate solutions for uni- 
axially or fully anisotropic d — 3 physical models. The global phase diagrams, with d — 2 and d — 1 
to d — 3 crossovers, are obtained for Ising, XY magnetic models and percolation systems, including 
crossovers from algebraic order to true long-range order. 
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I. INTRODUCTION 

Spatially anisotropic systems greatly enrich our expe- 
rience of collective phenomena, as exemplified by high-T c 
superconducting materials, in which the couplings along 
one direction are much weaker than those in the per- 
pendicular plane. Anisotropic systems are also intrigu- 
ing from a conceptual point of view, since vastly differ- 
ent critical phenomena are known to happen in differ- 
ent spatial dimensions, whereas, between d-dimcnsional 
systems stacked along a new direction, even the weak- 
est coupling, while not affecting the critical tempera- 
ture, induces (d + l)-dimensional critical behavior. Cal- 
culational results that yield the global phase diagram of 
anisotropic systems and thus provide a unified connected 
picture of the various anisotropic and isotropic behav- 
iors at different subdimensions and at the full dimen- 
sion have been rare and mostly confined to d = 2. In 
the present study, we obtain global phase diagrams for 
a variety of anisotropic d — 3 systems: Ising magnetic, 
XY magnetic, and percolation systems. Anisotropy along 
one direction (uniaxial) and full anisotropy, in which the 
couplings along each direction are different, are stud- 
ied, yielding global phase diagrams. We use hierarchi- 
cal models, which yield exact renormalization-group so- 
lutions [3, 0|. Thus, the construction of hierarchi- 
cal lattices that incorporate correct dimensional reduc- 
tions is an important step of the study. The exact 
solutions of hierarchical models can simultaneously be 
considered approximate position-space renormalization- 
roup solutions of models on naturally occurring lattices 
1]. The method developed in this study will be em- 
ployed to extend, from isotropic to anisotropic systems, 
the renormalization-group solutions of the tJ and Hub- 
bard models of electronic conduction 0, [H, Q ■ 



II. ANISOTROPIC HIERARCHICAL LATTICES 

Hierarchical lattices are constructed by repeatedly self- 
imbedding a graph. These provide exactly solvable 
models, with which complex problems can be studied 
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FIG. 1: The parallel, mutual imbeddings of these two graphs 
provide an anisotropic d = 2 hierarchical lattice. If either of 
the couplings (K x ,K y ) is set to zero, the remaining coupling 
constitutes a one- dimensional lattice. When the couplings are 
of equal strength, K x = K y , the two directions, represented 
by the two imbedding sequences, are equivalent and the lattice 
is isotropic d — 2. This lattice will be referred to as A^. 



and understood. For example, frustrated Q, spin-glass 
Q, random-bond @ and random- field [i"o| . Schrodinger 
equation [ll[ , lattice- vibration [l2j , dynamic sca ling [l3T | , 
aperiodic magnet complex phase diagram [l5| | and 
directed-path systems, etc., have been solved on hi- 
erarchical lattices. 

In this study, we construct anisotropic hierarchical 
lattices by the parallel, mutual imbedding of several 
graphs. In each imbedding step, b and b d respectively 
are the length and volume rescaling factors. We illus- 
trate the method by the simplest case of the anisotropic 
d = 2 lattice, before moving on to the uniaxially or 
fully anisotropic d = 3 lattices. The parallel, mutual 
imbeddings of the two graphs shown in Fig.l provide 
an anisotropic d = 2 hierarchical lattice. If either of 
the couplings (K x , K v ) is set to zero, the remaining cou- 
pling constitutes a one-dimensional lattice. When the 
couplings are of equal strength, K x = K y , the two direc- 
tions, represented by the two imbedding sequences, are 
equivalent and the lattice is isotropic d = 2. This lattice 




FIG. 2: The parallel, mutual imbeddings of these two graphs 
provide a uniaxially anisotropic d — 3 hierarchical lattice. If 
the coupling K z is set to zero, the coupling K xy constitutes 
an isotropic two-dimensional lattice. If the coupling K xy is 
set to zero, the coupling K z constitutes a one- dimensional 
lattice. When the couplings are of equal strength, K xy = K z , 
the z direction, represented by the first imbedding sequence, 
and the x,y directions, represented by the second imbedding 
sequence, are all equivalent and the lattice is isotropic d — 3. 
This lattice will be referred to as U3. 

will be referred to as A^. 

It is thus seen that generally our requirements in the 
construction of anisotropic hierarchical lattices are (1) 
the proper reduction to the lower dimension when one (or 
more, see below) of the couplings is set to zero and (2) 
the restitution of an isotropic lattice when the couplings 
are of equal strength. 

The parallel, mutual imbeddings of the two graphs 
shown in Fig. 2 provide a uniaxially anisotropic d = 3 
hierarchical lattice. If the coupling K z is set to zero, the 
coupling K xy constitutes an isotropic two-dimensional 
lattice. If the coupling K xy is set to zero, the coupling 
K z constitutes a one-dimensional lattice. When the cou- 
plings are of equal strength, K xy = K z , the z direc- 
tion, represented by the first imbedding sequence, and 
the x, y directions, represented by the second imbedding 
sequence, are all equivalent and the lattice is isotropic 
d = 3. This lattice will be referred to as U3. In Fig. 3, 
the parallel, mutual imbeddings of the top graph with 
either one of the following graphs also provide uniaxially 
anisotropic d = 3 hierarchical lattices. If the last graph 
is used, isotropy is not restored when K xy = K z ; this 
lattice is nevertheless included, for comparison, in our 
study. These lattices, differentiated by the choice of the 
second imbedding graph, will be respectively referred to 
as U 3a , U 3b , U 3c . 

A fully anisotropic d = 3 hierarchical lattice is pro- 
vided in Fig. 4 by each shown imbedding in parallel with 
the two imbeddings obtained by permuting K x (full line) , 
K y (dashed), and K z (dotted). If any one the couplings 
K u is set to zero, the remaining two couplings constitute 
an anisotropic two-dimensional lattice. If any two of the 




FIG. 3: The parallel, mutual imbeddings of the top graph 
with either one of the following graphs provide uniaxially 
anisotropic d — 3 hierarchical lattices. If the last graph is 
used, isotropy is not restored when K xy = K z . These lat- 
tices, differentiated by the choice of the second imbedding 
graph, will be respectively referred to as U3 a , Uzb, Us c - 

couplings are set to zero, the remaining coupling consti- 
tutes a one-dimensional lattice. When the couplings are 
of equal strength, K x — K y — K z , the three directions, 
represented by the three imbedding sequences, are equiv- 
alent and the lattice is isotropic d — 3. These lattices will 
be referred to as A 3a and A 3 b. 

The anisotropic systems that we study are located 
on the anisotropic lattices constructed above. These 
hierarchical models admit exact renormalization-group 
solutions, with recursion relations obtained by decima- 
tions in direction opposite to their construction direc- 
tion. The exact solutions of hierarchical models can si- 
multaneously be considered approximate position-space 
renormalization-group solutions of models on naturally 
occurring lattices. In fact, the recursion relations ob- 
tained for the models below correspond to Migdal- 
Kadanoff (l7l [Lsj approximate recursion relations, which 
are hereby generalized to anisotropic systems. 
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FIG. 4: Each imbedding shown in this figure, in parallel with 
the two imbeddings obtained by permuting K x (full line), 
K y (dashed), and K z (dotted), yields a fully anisotropic d — 3 
hierarchical lattice. If any one the couplings K u is set to zero, 
the remaining two couplings constitute an anisotropic two- 
dimensional lattice. If any two of the couplings are set to zero, 
the remaining coupling constitutes a one-dimensional lattice. 
When the couplings are of equal strength, K x = K y — K z , the 
three directions, represented by the three mutual imbedding 
sequences, are equivalent and the lattice is isotropic d — 3. 
These lattices will be respectively referred to as A^a and ^31,. 



III. ANISOTROPIC ISING MAGNETS 

The Ising model is defined by the Hamiltonian 
-0H 



S jS ; 



(1) 



where, at each lattice site i, Si = ±1, and < ij > u 
denotes summation over bonds of type u. The various 
decimations in the models are composed of two elemen- 
tary steps, K = K u + K v for bonds in parallel and 
K = tanh~ (tanni^ + tanhA'„) for bonds in series, 
where K is the effective coupling of the combined bonds. 

The phase boundaries for the Ising model on the d = 2 
anisotropic hierarchical lattices A 2 ,A 3a , and A 3b (set- 
ting K z — in the latter two) are given in Fig. 5, along 
with the exact result for the anisotropic square lattice 
[l9l ]. The renormalization-group flows are indicated on 
the phase boundary of the hierarchical models. The fixed 
point occurs at isotropy, K x = K y , to which the d = 2 
anisotropic critical points flow, thereby sharing the same 
critical exponents. 

The phase boundaries for the Ising model on the 
d = 3 uniaxially anisotropic hierarchical lattices 
U3, U 3a , U 3b , U 3c ,A 3a , and A 3b (setting K x = K y ) are 
given in Fig. 6. The exact phase transition points for the 
square and cubic [2(| lattices are also shown. For 
each model, the phase transitions at d = 1 (at infinite 
coupling) and d — 2 cross over to d = 3 criticality, which 
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FIG. 5: The phase boundaries for the Ising model on the d — 
2 anisotropic hierarchical lattices A £ ,Ao, a , and A^b (setting 
K z = 0) and the exact result for the anisotropic square 
lattice (s). The arrows and stars respectively indicate the 
renormalization-group flows and fixed points. 



is thus universal for all d = 3 anisotropic and the d = 3 
isotropic cases. 

The phase boundary surface for the Ising model on the 
d = 3 fully anisotropic hierarchical lattice A 3 b is given in 
Fig. 7. The dashed lines on the planes are the exact d = 2 
solutions for the square lattice fl9| . Again, all points on 
the critical surface of the d — 3 fully anisotropic model 
flow onto the fixed point located at isotropy, thereby 
sharing its critical exponents. The critical exponents 
found for this model are j/t = 0.69, yn = 1-68 for d = 2 
(for the square lattice y T — \ 1 y H = 1.875 [19 1 ]) and 
yr = 0.92, yn = 2.20 for d = 3 (for the cubic lattice 
Kr = 1.59, j/jf = 2.50 0,11). 
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FIG. 6: Phase boundaries for the Ising model on 
the d = 3 uniaxially anisotropic hierarchical lattices 
U3,U: ia , U-ib, U-jc, A 3a , and A 3b (setting K x = K y ). The exact 
phase transition points for the square JTtJ and cubic [2(| lat- 
tices are shown by the black circles. For each model, the phase 
transitions at d — 1 (at infinite coupling) and d = 2 cross over 
(as shown for As a ) to d = 3 criticality, which is thus universal 
for all d = 3 anisotropic and the d = 3 isotropic cases. The 
d = 2 fixed point of A^b is not marked by a star, since it co- 
incides with the square lattice exact transition point, marked 
by the black circle on the horizontal axis. 



IV. ANISOTROPIC XY MAGNETS 

The XY model is denned by the Hamiltonian 

-PH = ^ .4 ■ S, = ^ J„ ^ 008(0; - 6 S ), (2) 

u (ij)u u (ij)u 

where at each lattice site i, S; is a unit vector con- 
fined to the xy plane at angle 9i to the x axis and 
< ij > u denotes summation over bonds of type u. Un- 
der renormalization-group transformations, the coupling 
between nearest-neighbor sites takes the general form of 
a function V u {9i — 9j). The various decimations in the 
models are composed of two elementary steps, 

V =V U + V v and (3) 

Vie, -6 k )= In / d6j exp[K (fc - 6 3 ) + V v [9, - d k )}, 
Jo 

respectively for bonds in parallel and in series, where V 
is the effective coupling of the combined bonds. In terms 








FIG. 7: Phase boundary surface for the Ising model on the 
d — 3 fully anisotropic hierarchical lattice ^36. The dashed 
lines on the planes are the exact d — 2 solutions for the square 
lattice [3. 

of Fourier components, 

U(s) = f V ^ se exp[V u (6) - K(0)], (4) 
Jo 2?r 

exp[K(#)-K(0)] = J2e- ise fu(s), 

s 

Eqs.(3) respectively are 

f(s) =^2fu(p)fv(s-p) and (5) 
p 

f(s) =U(s)f v (s), 

in a form that is more conveniently followed in our cal- 
culations. The phase boundaries for the XY model on 
the d = 3 uniaxially anisotropic hierarchical lattices 
U3, U^a, U3b, and Us c are given in Fig. 8. The exact phase 
transition points for the square [23| and cubic [24| lattices 
are also shown. 

In d — 2, namely along the horizontal axis, above a 
critical interaction strength marked by the squares on 
the figure, the systems exhibit algebraic order [25|, [2(| 
\27\ : The starting Hamiltonian [Eq.(2)] flows to a Villain 
potential [H, 

f v (s) =Aexp(- s 2 /2Jy), (6) 

exhibiting a fixed-line behavior parametrized by Jy. This 
corresponds to a system without true long-range order, 
namely with zero magnetization, but infinite correlation 
length and algebraic order in which the correlations de- 
cay as a power law. In d = 3, true long-range order oc- 
curs: points in the ferromagnetic phase renormalize to a 
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FIG. 8: Phase boundaries for the XY model on the d=3 uni- 
axially anisotropic hierarchical lattices Us,Usa,U3h, and Use 
The exact phase transition points for the square [2^] and 
cubic [2J] lattices are shown by the black circles. For each 
model, the phase transitions at d = 1 (at infinite coupling) 
and d = 2 (onset of algebraic order) cross over to d = 3 crit- 
icality, which is thus universal for all d = 3 anisotropic and 
the d = 3 isotropic cases. The onsets of effective algebraic 
order in d = 2 are marked, for models U3, U3 a with the open 
square and for models U3b,Us c with the full square. 



delta function potential; points on the phase boundaries 
renormalize to single true fixed potential, shown in Fig. 9, 
differing from the Villain potential as also seen on the fig- 
ure. The behavior here for d = 2 is not true fixed-line 
behavior. After tens of thousands of renormalization- 
group iterations (corresponding to a scale change factor 
of 210,000), the Villain potential decays [27j to a disor- 
dered sink with V(9) max — V(6) m i n < 10~ 4 . The sharp 
change in the necessary number of iterations, as seen in 
Fig. 10, indicates the onset of effective algebraic order. 

As seen in Fig. 8, for each XY model, the phase tran- 
sitions at d — 1 (at infinite coupling) and d — 2 (onset 
of algebraic order) cross over to d — 3 criticality, which 
is thus universal for all d = 3 anisotropic and the d = 3 
isotropic cases. 



V. ANISOTROPIC PERCOLATION 

Anisotropic percolation is defined such that on each 
connection of direction u, a bond exist with probability 
p u . The various decimations in the models are composed 
of two elementary steps, p — p u Pv+Pu(l—Pv)+Pv({—Pu) 
for connections in parallel and p = p u p v for connections 
in series, where p is the effective connectedness probabil- 
ity of the combined connections. The phase diagram for 
percolation on the d = 2 anisotropic hierarchical lattice 




j 1 i_ 

45 90 135 180 
6 (deg) 



FIG. 9: Fixed potentials attracting the critical surface of the 
d — 3 XY model. For comparison, appropriately normalized 
Villain potentials are also shown. 




1/J Villain 

FIG. 10: Number of renormalization-group iterations nec- 
essary, in the d = 2 XY model, for the Villain poten- 
tial parametrized by Jv to decay to a disordered sink with 
< 10~ 4 . The squares indicate the onset 
of effective algebraic order. 



(A2) is given in Fig.ll. The percolation fixed point oc- 
curs at isotropy, p x — p y , to which the d — 2 anisotropic 
percolation onsets flow, thereby sharing the same crit- 
ical exponents. The phase boundaries for percolation 
on the d = 3 uniaxially anisotropic hierarchical lattices 
U3,U3a,U3b, and U^c are given in Fig. 12. The perco- 



6 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



non 




percolating ^ 



6* 



0.0 



0.2 



0.4 



0.6 0.8 



1.0 



FIG. 11: The phase boundaries for percolation on the d = 2 
anisotropic hierarchical lattice (A2). The renormalization- 
group flows are indicated on the phase boundary of the hier- 
archical model. The onset of percolation for the square lattice 
[29l | is shown by the black circle. 



lation poin ts for the cubic and square lattices are also 
shown 29]. For each model, percolation onset at d = 1 
(at p z = 1) and d = 2 cross over to d — 3 percolation 
onset, which is thus universal for all d = 3 anisotropic 
and the d = 3 isotropic cases. 
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FIG. 12: Phase boundaries for percolation on the d = 3 uni- 
axially anisotropic hierarchical lattices Uz,Uza,U3b, and Us c . 
The onsets of percolation for the square and cubic lattices 
[29j are shown by the black circles. For each model, the per- 
colation onset at d — 1 (at p z = 1) and d — 2 cross over to 
d = 3 percolation onset, which is thus universal for all d = 3 
anisotropic and the d = 3 isotropic cases. 
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